Preprint typeset in JHEP style. - HYPER VERSION 



DCPT-02/39 



A new analytic approach to physical 
observables in QCD 



Wajdi Gaddah 

Centre for Particle Theory 
Department of Mathematical Sciences 
University of Durham 
South Road 

Durham, DH1 3LE, England 



E-mail: H . A . GaddahOdurham .ac.uk 



Abstract: An analytic ghost-free model for the QCD running coupling a(Q 2 ) is 
proposed. It is constructed from a more general approach we developed particularly 
for investigating physical observables of the type F(Q 2 ) in regions that are inaccessi- 
ble to perturbative methods of quantum field theory. This approach directly links the 
infrared (IR) and the ultraviolet (UV) regions together under the causal analyticity re- 
quirement in the complex Q 2 — plane. Due to the inclusion of crucial non-perturbative 
effects, the running coupling in our model not only excludes unphysical singularities 
but also freezes to a finite value at the IR limit Q 2 = 0. This makes it consistent 
with a popular phenomenological hypothesis, namely the IR freezing phenomenon. 
Applying this model to compute the Gluon condensate, we obtain a result that is in 
good agreement with the most recent phenomenological estimate. Having calculated 
the 13— function corresponding to our QCD coupling constant, we find that it behaves 
qualitatively like its perturbative counterpart, when calculated beyond the leading 
order and with a number of quark flavours allowing for the occurrence of IR fixed 
points. 
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1. Introduction 

A fundamental problem in the theory of particle physics lies in the description of 
hadron interactions in the infrared region. The property of asymptotic freedom 
|]] in QCD allows us to investigate the interactions of quarks and gluons at short 
distances using the standard perturbation theory. However, there is a number of 
phenomena whose description is intractable in perturbation theory, for example quark 
confinement and gluon and quark condensates. Hence, the search for calculational 
techniques that go beyond conventional perturbation theory remains essential (not 
only in QCD but in other field theories as well). 

In standard perturbative QCD, for any number of quark flavours n/, the running 
coupling constant (to leading order): 

a ^ Q2)= j B M^m " ith * =n -&>- (L1) 

diverges at a small mass scale Q = A, creating the so-called Landau ghost-pole prob- 
lem. Taking next loop corrections into account does not alter the essence, and leads 
only to additional branch cuts along the positive real axis Re{Q 2 } > 0. This prob- 
lem prevents the use of perturbative expansion at small momentum transfers Q ~ A 
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and, in addition, generates infrared (IR) renormalon singularities on the positive real 
axis of the Borel parameter, destroying attempts to sum up the perturbative series 
[0, |3|, £|, |j. Generally speaking, as Q 2 comes below or near A 2 , non-perturbative 
effects become the most dominant and the perturbative expansion becomes useless. 

Another indication that the perturbative formalism is incomplete, and cannot 
describe the low energy physics unless it is supplemented by non-perturbative cor- 
rections, comes from considering its analyticity structure in the complex Q 2 -plane 
0; the upshot being that any QCD observable, which depends on a spacelike 1 
momentum variable Q 2 , is expected to be an analytic function of Q 2 in the entire 
complex plane except the negative real (timelike) axis. Singularities on the timelike 
axis are meaningful since they correspond to production of on-shell particles, while 
their existence on the spacelike axis is non-physical as this would violate causality 
(i.e. causal analyticity structure) [7|. For example, if a generic QCD observ- 
able were calculated to leading order, it would depend on ap^(Q 2 ), inheriting the 
Landau- singularity at Q 2 = A 2 . This singularity is obviously non-physical as it lies 
on the positive real axis Ke{Q 2 } > 0. From this point of view, causality constrains 
the Q 2 — dependence of physical observables, and is consistent with perturbative re- 
sults only if the coupling constant is non-singular on the entire Q 2 -plane with the 
exception of the timelike axis Re{Q 2 } < 0. 

A small number of models for the QCD running coupling a(Q 2 ) have been pro- 
posed [|, 0, ||, [|, [10] to comply with the causality condition, i.e. the requirement of 
analyticity in Q 2 . However, the issue of which of these is the most realistic is still a 
moot point. Quite recently, a good attempt has been made, by Shirkov and Solovtsov 
0, to devise a model for the running coupling that is completely free of singularities 
in the IR region. In this approach, the analytic running coupling is defined via the 
so-called Kallen-Lehmann spectral representation as: 



where aipij, is the coupling constant obtained from perturbation theory (PT) at the 
n-loop approximation. To leading order, the corresponding spectral density reads as: 



1 In this paper, the metric with signature (-1,1,1,1) is used so that Q 2 > corresponds to 
spacelike (Euclidean) 4-momentum transfer. 





with the n-loop spectral density (x) given by 
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Inserting this into ( 



OP gives the one- loop analytic spacelike coupling constant: 



a 



(i) 



(Q 5 



47T 
To 



_ln(Q 2 /A 2 ) 1 -Q 2 /A 2 _ 



1.5) 



This expression is consistent with the causality condition. The first term on the 
right-hand side of ( |1.5|) preserves the standard UV behaviour whereas the second term 
compensates for the ghost-pole at Q 2 = A 2 . Obviously the second term in (|1.5|) comes 
from the spectral representation and enforces the required analytic properties. Thus, 
it is essentially non-perturbative. However, as mentioned in Ref.pj, as this term 
introduces a l/Q 2 correction to ap^(Q 2 ) at large Q 2 it would make ( |1.5| ) unsuitable 
as an input quantity for observables that are proportional to the running coupling 
(in leading order) but are not expected to have l/Q 2 corrections (such as the e + e~ 
total cross section) because in such cases the unwanted l/Q 2 correction would have 
to be artificially cancelled. By construction, ( |T3|) does not include any adjustable 
parameters other than the QCD characteristic mass scale A. Thus, one might prefer 
to have a model with free extra parameters that would allow for the fitting of the 
coupling to lower energy experimental data and hence further improve perturbative 
results at low Q 2 . An example of a model of this type is suggested in Ref. 0, where 
the one-loop coupling constant assumes the form: 



47T 



* (1) (Q 2 ) = f 



Q 2p /A 2p 



[Q 2 p/A 2 p + C p 



P 



\n(Q 2 P/A 2 P + C p ) 



1.6) 



with C p > 1. 



For Q 2 > 



this expression recovers the perturbative asymptotic 



form ( |1 . 1|) . Although (L6) is a ghost-pole free expression, being analytic for all 



Q 2 e [0, oo) , it contains unphysical singularities on the complex Q 2 — plane (for all 
p > 1) and hence does not comply with causality. A careful analytic study of non- 
perturbative effects is always important in QCD because it provides more reliable 
and useful information about the IR region, which is inaccessible to perturbative 
methods. 

The main purpose of this paper is to construct a new analytic approach based 
on a different and more general methodology to that of Shirkov and Solovtsov, which 
improves perturbative results outside the asymptotic domain and respects the causal 
analyticity structure in, at least, the right half of the complex Q 2 — plane. Having 
achieved this, we apply our method to the QCD coupling constant to solve the 
Landau ghost-pole problem without altering the correct perturbative behaviour in 
the UV region. Then, we show that the resultant QCD coupling in this way freezes 
to a finite value at the origin (i.e. Q 2 = 0), supporting the IR freezing idea which has 
long been a popular and successful phenomenological hypothesis ||TT| , |T~2j] . Moreover, 
from the viewpoint of the new background field formalism |3|, Q our approach provides 
better estimates for the IR fixed points than that of Shirkov and Solovtso. 
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The layout of this paper is as follows: In section 2, we discuss our analyticization 
procedure on general grounds, proposing a new integral representation for physical 
observables that depend on a spacelike momentum variable (squared) Q 2 . In the 
next section, we present a full derivation of our analytic coupling constant, giving 
a new expression that depends on a free extra parameters A e which can be used to 
tune the coupling to lower energy experimental data. An approximation scheme for 
estimating the value of A e is included in section 4. In section 5, we carry out a com- 
prehensive comparison between the predictions in our approach and those estimated 
by other theoretical methods, which include conventional perturbation theory, op- 
timized perturbation theory, background field formalism and the analytic approach 
of Shirkov and Solovtsov. In this section, we test our model on a fit-invariant IR 
characteristic integral extracted from jet physics data. For further applications, we 
use our analytic coupling constant in the instanton density derived from the dilute 
instanton-gas approximation to estimate the gluon condensate. Then, we show that 
the result obtained is in good agreement with the value phenomenologically estimated 
from the QCD sum rules. Finally, we end this section by calculating the (3— function 
which corresponds to our analytic coupling constant and compare its behaviour with 
the perturbative counterpart. We further show that the /3— function in our approach 
and that in higher-loop order perturbation theory behave qualitatively the same in 
the range of rif— values that allows perturbation theory to obtain IR fixed points. 
The last section is devoted to our conclusions. 



2. Analyticization Procedure 

We will illustrate our method with the following pedagogical example. Consider 
a physical observable F depending on a positive variable q with the dimensions of 
energy squared. Then, it follows from the principle of causality, mentioned earlier, 
that F(q) can be analytically continued to the complex plane excluding the negative 
real axis. As we shall show, this will allow us to reconstruct F(q) from its high energy 
behaviour by using the contour integral representation: 

F(?,A) = — / F(k) dk , (2.1) 

2m Jc uv k-q 

where Cuv is a very large incomplete circle, in the complex k— plane, with radius 
Ruv an d centre at the origin, beginning just below the negative real axis and ending 
just above and q & {Re k : < Re k < Ruv} as depicted in Fig.(|l]). 

We will show that F(q) is recovered as the limit as A — > oo of F(q,\), and 
this depends mainly on the ultraviolet behaviour of F(q) when we take Ruv to be 
sufficiently large. In order to achieve this goal, we exploit the analyticity of F(k) in 
the complex fc-plane, using Cauchy's theorem to replace the UV-boundary Cuv m 
( |2.1| ) with a small circle Cq around q and a keyhole shaped contour Ck, depicted in 
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Figure 1: Sketch of the integration contour Cjjy i n the cut complex k— plane. The thick 
line represents a branch cut and the cross indicates the location of q. 
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Figure 2: Sketch of the keyhole shaped contour Ck- The thick line represents the branch 
cut Re k < 0. Here, e is an innnitesimally small radius and 9 S (— tt, 7r] . 



Fig.(0), surrounding the cut and not connected to Cq. For this deformed contour, 
Cauchy's integral formula over the circle Co yields F(q) which together with the 
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integral over the keyhole contour Ck transforms (|2.1|) into the form: 



F(q,\)= F(q) + I Ck , 



[2.2) 



where 



l c k 



1 

2vri 



F(k) dk. 



(2.3) 



<c k K-q 

If we restrict attention to those functions F(k) which are finite on the line of discon- 
tinuity along the negative real /c-axis, e.g. F(k) oc \f~k~, then we can write Ic k more 
explicitly as: 



l c k 



Ruv g— A(r+g) 



r + q 



A(r) dr + 



e c 



-M 



2tt 



with 



and 



A(r) 



W,q) = - 



2m 



[Fire 1 *) - F(re 



e Aecos(0) f e (0, q) d6, (2.4) 



(2.5) 



e i\esm(d) P^qM} e -iAesin(6») p^Q~ ie ^ 



e — q e 



-w 



e — q e 



i6 



(2.6) 



Owing to the discontinuity of F(A;) across the cut, given by 2ni A(|fc|) , the sum of 
the line integrals on each side of the branch cut does not cancel. Since the region of 
analyticity of F(k): 



{k = re ie : r > , -tt < 9 < n}, 



(2.7) 



includes a portion of the real axis on which F(k) is real, it follows that F(k) satisfies 
the Schwarz Reflection Principle ||14|| : 



F*(k) = F{k*). 

This allows us to rewrite A(r) and f e (6 l ) in the more compact forms: 



AM 



TT 



Im{ F{re i7T ) } , 



and 



U9,q) = Re 



-jiAe sin(0) 



F(ee 



i0\ 



q e 



(2.{ 



(2.9) 



(2.10) 



which imply that is a purely real quantity as it should be. 

In the following, we will investigate the behaviour of Ic k at sufficiently large 
values of A. Since F(re z9 ) is assumed finite throughout the whole region S an , we 
take its magnitude to be bounded on the infinitesimal circle \k\ = e, see Fig.(|2|), such 
that: 

\F(ee id )\<M, for - tt < 9 < n , (2.11) 
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where M is a positive constant which may depend on e. Using this to estimate the 
upper bound of the second term in ( |2.4j ), we find that: 



6 c 



-Xq 



2vr 



e Xecos(e) ^ 



< 
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(2.12) 



If we consider only those functions F for which M can be chosen such that lim eM = 

then the second term in (2~ 3) will vanish in this limit and accordingly Ic k reduces to 
the form: 



ic k (q, A) 



#E/V g-A(r+g) 



r + g 



A(r) dr, 



g > 



(2.13) 



which is more suitable for further calculations. For the case under consideration, 
since A(r) has an upper bound N > 0, i.e. |A(r)| < N , in the interval < r < Ruv, 
we deduce that: 



-Xq 



< N 



Xq 



(1 



-XRuv 



)■ 



(2.14) 



So for large values of A, the function Ic k (q, A) decays rapidly to zero, and: 

lim I Ck {q,X) = 0. (2.15) 

A— >oo 

In general, ( |2.15| ) remains true even if the finiteness restriction on F(k) along the 
negative real fc-axis, together with the constraint limeM = 0, is not fulfilled. For 
example, if F(k) has a finite number of singularities on the negative real axis then 
we can avoid these points by keeping the two edges of the keyhole contour Ck a small 
distance 5 away from the cut which they surround. In this way, we can rewrite ( |2.3| ) 
as: 



l c k 



Ruv 



-X(x+q) 



g(x, g, A) dx + 



€ c 



-Ay 



2tt 



TT—r] 



TT+T] 



where e' = \/ e 2 — 5 2 with e > 5 , r\ = arctan (S/e') and: 



1 f e iXS ) 

g(x, q, A) = - Im \ — F(-x +iS) \ . 

7i I x + q — to J 



x + q — iS 

Since F is bounded on the contour CV , we have: 



< 



N' 



-Xq 



7T \(q + £' -8) I 



g— Ac' q-^Rjjv 



+ 



eM{<K-rj) e'^i-") 



7T 



(2.17) 



(2.18) 



where N' is the upper bound of F(—x + i5) for all x G [e', Ruv] ■ As both N' and M 
are A-independent, in the limit as A — > oo we obtain ( p,15| ) as claimed above. So, by 
taking the limit as A — > oo in ( |2.2j ) we deduce from ( |2.15| ) that: 

lim F(q, A) = F(q) . (2.19) 
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To have some notion of the functional dependence of F(q, A) on A, for a fixed 
q > 0, consider the simple example where the imaginary part of F(k) assumes a 
constant value /i along the upper edge of the cut, e.g. F(k) oc In A;. In this case, we 
have A(r) = /U /7r and a straightforward integration of ( |2.13| ) yields: 



l c k (q, A) 



/'o 



7T 



E 1 (Xq)-E 1 (X(q + R uv )) 



(2.20) 



where E\{x) denotes the exponential integral function [T5| defined by: 



Et(x) 



oo -t 

— dt 
t 



"7 



ln(x) - 



n=l 



-1) 



ra! n 



x 



for x > 



(2.21) 



here 7 = 0.577215 is the well-known Euler's constant. By inserting ( |2.20| ) into ([27 
and sending Ruv — > 00, we immediately arrive at an explicit solution: 



F(q,X) = F(q)+^E 1 (Xq). 

7T 



(2.22) 



As illustrated by the graph in Fig.(||), the variation of F(q,X) with respect to A 
depends largely on both the sign of fio and the behaviour of Ei(Xq). With increasing 
A, we see that F(q, A) evolves from —00, for fi < 0, or +00, for fi > 0, in a 
continuous progression until it settles down to the value F(q) at large values of A. 



Moreover, we discover from (|2.22|) that F(q, A) converges to F(q) faster as q increases. 



To summarise, by exploiting the analyticity of F(k) in ( |2.1| ), we can obtain F(q) 
as lim^oo F(q, A), and if we take the radius Ruv of the contour Cuv to be sufficiently 
large the integrand depends on the large q behaviour of F. So far, our consideration 
has been restricted to the case in which F(q), as a true function for the observable 
in question, complies with the principle of causality. In practice, however, we do not 
know the function F(q) exactly, even for large q, but for a theory with asymptotic 
freedom we can approximate the large q behaviour of F(q) using perturbation theory 
for q in some high energy domain Dp? so that Fpx(q) — F(q) for q G Dp?- Thus, 
for large enough Ruv, we can safely use F PT {k) as a reasonable substitute for F{k) 
in our contour integral ( [2.1Q , where now k is any point on Guv- This gives the 
approximation Fpr(q,X) ~ F(q,X), but taking A to infinity in this would simply 
reproduce the perturbative estimate for F(q), leading to nothing new. Now Fpx{q), 
as a perturbatively calculated quantity, often suffers from a serious defect by having 
singularities in the complex plane away from the negative real axis. Good examples 
of this are the effective coupling constant and the running mass in perturbative QCD. 
These extra singularities are unphysical and do not exist in the true function F(q). 
Hence, they should be eliminated in our approximation in order to reinstate causality. 
So, if instead of using the limiting value of F PT (q, A) we were to use a finite value of 
A, say Aq, then we might have an expression in which these unphysical singularities 



S 



F(q, A,) iL , 



F(q) 








Figure 3: The variation of .F(<7, A) versus A at a fixed value q . The lower and the upper 
curves describe the cases in which no < and > 0, respectively. The dashed line 
represents F(q). 



are absent and we will give an example in which this happens. Furthermore, taking a 
finite value of A in Fpr(q, A) can give an improvement over the perturbative estimate 
in favourable circumstances for the following reason: Viewed as a function of A with 
q held fixed, F PT (q, A) is a better estimate of the true function F(q, A) for small A 
than it is for large A since if we scale the integration variable in fl2.1| ) then: 

F PT (q,\) = — / ^_^j?pr(VA)dfc. (2.23) 

2™ J\ k \=\R uv [k - Xq) 



For large enough Ruy we can replace the contour of integration in ( |2.23| ) hy \k\ = R, 



Wv 
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since the difference this makes is a small contribution from the negative real axis: 

" R uv p— (k+Xq) ( 

Im{ F PT (ke i7T /X) } dk 



jo — - 

7T 



A Rjj\ 



(k + Xq) 



2.24) 



which is exponentially damped. Now we can see that the smaller the value of A is the 
larger is the argument of Fpt in the integrand of (|2.23|) , and for an asymptotically 
free theory this gives a better approximation to F(q,X). So, FpT(q,X) interpolates 
between the true quantity F(q, X) at small A and the perturbative estimate F PT (q) 
at large A. Under favourable circumstances A can be chosen sufficiently small that 
F PT (q, A) provides a good enough approximation to the exact function F(q, A), and 
sufficiently large that we are estimating the large A behaviour of the latter, i.e. F(q). 

We will now discuss how to choose an appropriate value A = Ao for FpT(q,X), 
depending on the behaviour of Fpt, and pinpoint the most favourable situations in 
which F PT (q, A ) gives a better estimate for F(q) than F PT (q) does. An illustrative 
example is perhaps the best way to convey these ideas. Consider the simple case in 
which Fpx{q) suffers only from one infrared singularity at q = A 2 of the simple pole 
type: 

u(q) 

q- A 2 

where U(q) is an analytic function of q in the entire complex q— plane excluding the 
negative real axis (e.g. the principle value of log). Then, define the UV region where 
the perturbation theory of this toy model is trustworthy by: 

D PT = {q ■■ q > A 2 } . (2.26) 

Replacing F and F with Fpt and Fpt-, respectively, in our integral formula ( |2.1| ) 
and then integrating in the usual way, we obtain: 

F PT (q, A) = F PT (q) - e^"^ + H(q, A) for q ? A 2 , (2.27) 

q — A z 



Fp T (q) = i^- 2 , (2.25) 



and 
where 

with 



F PT (A 2 , A) = A U(A 2 ) + U'(A 2 ) + H(A 2 , A) for q = A 2 , (2.28) 

[•Ruv p-A(r+g) 

H(q,X)= ■ A PT (r)dr, q>0, (2.29) 

Jo r + Q 



1 Im{ U(r e ln ) } 
7T r + A 2 

Here, we have Ruv ^ A 2 . In this model, (|2.27f) allows us to modify the perturbative 



M')'-; '.! , 1 ' . (2-30) 



expression FpT(q) through the variable A. Taking A to infinity in fl2.27]) would only 
reproduce F PT (q)- 

lim FpT(q, A) = Fpp(q) , (2.31) 

A— >oo 
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on the other hand, keeping A fixed at a finite value A would change the behaviour 
of FpT{q), especially in the vicinity of A 2 , by an infrared correction term: 

U(A 2 ) 



Ym(q, A c 



3 -A (g-A 2 ) 



+ H(q, A c 



(2.32) 



q- A 2 

which would, in turn, remove the IR singularity in F PT (q). So, to reinstate the causal 
analyticity structure violated by the presence of the simple pole singularity on the 
positive real axis, we would simply retain Tm(q, Ao) in our formulation. For a finite 
value of A, the continuity of F PT (q, A) at q = A 2 : 



lim F PT (q,X) = F PT (A 2 , A) 

q^A 2 



(2.33) 



follows directly from (g]27]) and (|2~28|) . 

Now, we shall study this example further by investigating the possibility of find- 
ing a proper value A = A that would allow FpT^q, Ao) to match the exact value of 
the observable F(q) at the point q in question. At the same time, we shall consider 
the case in which this is not possible and show how to make the best choice of Ao 
that can improve perturbative predictions. To explore straightforwardly the crite- 
rion of selecting Ao, at a fixed value of q, we shall simplify our example further by 
assuming that Im {U(re l7T )} = u, with uo being some real constant. In this case, a 
straightforward calculation of ( |2.29p yields: 

UJ 



H(q, A) 



n(q - A 2 ) 



-A(q— A 2 ) 



E^A 1 ) - E X {X{A 2 + R uv )) 



E x {Xq) + E 1 (X(q + R uv )) } for q + A 2 



(2.34) 



and 



H(A 2 , A) 



uj 



71 



-AA 2 



-XRjjv 



- Ae 



AA 2 



A 2 A 2 + R uv 

- £ X (A(A 2 + R uv )) 



Ei(XA 2 



for q = A^ 



(2.35) 



Note that in the limit as q — > A 2 equation Q2.34 ) tends to ( p.35|) . Substituting ( p.34 ) 
into ( p.27| ) and then letting Ruv °o, w e obtain: 



F PT (q,X) = Fpx(q) - U(A 2 ) H Ei(\A ) - - E^Xq) e 

L 7T 7T 

which together with the first and the second derivatives: 



,A(<?-A 2 ) 



,-A(g-A 2 ) 
q- A 2 



Tx FMqA) 



U^ + ^E^XA 2 

71 



and 



UJ 



dX 2 ' 



71 



71 . 



F PT (q, A) = < (q - A 2 ) E X (XA ) + -U(A 



-\(q-A 2 ) 



-AA 2 



+ 



A 



-A(9-A 2 ) 



(2.36) 



(2.37) 



(2.38) 
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allows us to deduce and plot the graphs of the possible cases governing the variation 
of FpT^q, A) against A. These can be classified in the region q > A 2 according to the 
signs of the constants uj and U(A 2 ) in the following way: 

1. if uj > and U(A 2 ) > then: 

d d 2 

—F PT (q, A) > and -^Fpt(Q> X ) < > for A > , 
dX dX z 

which imply that with increasing A the curve of Fpj<(q, A) increases towards 
FpT^q) with a downward concavity as depicted in Fig.(£|.i.a) and Fig.(§.l.b); 

2. if uj > and U(A 2 ) < then there is a local maximum at A = X max , which can 
be deduced from: 



d_^ 
dX 

d 2 ^ 



^F PT (q,X) _ =0 =S- E 1 {\ max A 2 ) = -\U{A 2 )\ and 



dX- 



F PT (q, A) 



< 

A — ^m.n.T, 



Here, the way in which Fpp(q, A) approaches FpT^q) is different because Fpp(q, A) 
rises first to the maximum F PT (q, X max ) to fall next towards F PT (q) as A in- 
creases beyond X max as shown in Fig.(|.2.a) and Fig.(^.2.b); 

3. if uj < and U(A 2 ) < then: 

^F PT (q, A) < and ^F PT (q, A) > , for A > , 
dX dX z 

which indicates clearly that with increasing A the curve of FpT^q, A) decreases 
towards F PT (q) with an upward concavity as depicted in Fig.(^|.3.a) and Fig.(|5].3.b); 

4. if uj < and U(A 2 ) > then there is a local minimum at A = A m j n , which can 
be deduced from: 

d F PT (q,X) _ = £i(A mm A 2 ) = ^-U(A 2 ) and 



dX 
£_ 
dX? 



F PT (q, A) 



> 



Here, Fpr(q, A) decreases in the direction of the minimum FpT(q, A m j„) as A 
approaches X min from the left and begins to rise towards F PT (q) immediately 
after A passes X min as shown in Fig.(|5|.4.a) and Fig.([|.4.b). 

For each of these four cases, there corresponds two possibilities either 

(a) F{q)<F PT {q) or, 

(b) F(q) > F PT (q) . 
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(La) 




(Lb) 




(2.a) 




(2.b) 



Figure 4: The variation of FpT(q,\) and the expected F(q,\) versus A at a fixed value 
q > A 2 . In (l.a) and (2. a) the upper and the lower curves describe Fpp(q, A) and F(q, A) re- 
spectively. In (l.b) and (2.b) the upper and the lower curves describe F(q, A) and Fpx(q, A) 
respectively. 



Accordingly, we shall denote the case associated with either possibility by case(i.a) 
or case(i.b) with i — 1, 2, 3 or 4. These cases are illustrated in Fig.(|4]) and Fig.([5|). In 
these figures, we include for each case an illustrative hypothetical graph for F(q, A) 
based on our knowledge of the small and large behaviour of F(q, A), i.e. F(q, A) = 
Fpr{q, A) for sufficiently small A and F(q, A) = F(q) for large enough A. In applying 
our method we will need to know a priori whether Fpp{q) is less than or greater than 
F(q), either by input from an experiment or by a knowledge of the signs of higher 
order corrections to the perturbation theory in a regime where this can be trusted. 
Our choice of A for the different cases, at a fixed value q > A 2 , is as follows: 

1. in case(l.a), the ideal value of Aq is the one at which FpT{q, A) intersects with 
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F PT (q) 
F(q) 


v - 




\ 

> 


\ 


















(3.a) 




(3.b) 


F PT (q) 
F(q) 


\ 


f (tf ) 
F PT (q) 


\ 
























X q = A. m/'/i A. 

(4.a) 




A q A, min X 

(4.b) 



Figure 5: The variation of FpT(q,\) and the expected F(q,\) versus A at a fixed value 
q > A 2 . In (3. a) and (4.a) the upper and the lower curves describe Fpp(q, A) and F(q, A) re- 
spectively. In (3.b) and (4.b) the upper and the lower curves describe F(q, A) and Fpx(q, A) 
respectively. 



F(q) as shown in Fig.(|Jl.a); such a point is not precisely known although it 
may be estimated from a knowledge of higher order corrections to Fpt if these 
are available. We expect it to occur near the point of greatest curvature of 
F PT (q, A) as illustrated by the shaded region in Fig.([|.l.a). So for this case we 
might alternatively take A to be the value of A where the curvature is greatest. 

2. in case(l.b), the monotonicity of FpT{q, A) means that sampling this curve at 
any finite value of A will fail to improve the perturbative estimate. However the 
use of a large, but finite value will depart little from the perturbative estimate 
(as the asymptotic behaviour is approached rapidly) and will have the merit 
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of reinstating causality. So our best choice for Ao is the smallest A— value at 
which Fpr{q, A) is acceptably close to Fpr(q)- 

3. in case(2.a), the ideal value of Ao is the one at which Fpj<(q, A) intersects with 
F(q) as shown in Fig. ([|.2. a); such a point is expected to lie within a narrow 
domain, on the left of \ m ax, corresponding to a small arc segment just below 
Fpr{q) as illustrated by the shaded region in Fig.(|.2.a); such a point might be 
estimated from a knowledge of higher order perturbative corrections; however 
without further information it is not possible to see how to identify this region 
from a knowledge of the perturbative function FpT^q, A) alone. 

4. in case(2.b), the best choice for A is \ ma x, as illustrated in Fig.(|.2.b), and 
this leads to an improvement over the perturbative estimate. 

The cases depicted in Fig.(|5]) differ essentially by just a sign change, so a similar 
discussion of the appropriate choice of Ao can be made. 

To summarise, from a prior knowledge of whether the perturbative estimate 
overshoots or undershoots the true value of the observable F(q), and by plotting the 
shape of Fpx(q, A) against A we can ascertain which of the various possible cases 
occurs in a given situation. For cases (l.a), (2. a), (3.b) and (4.b), there exists a 
unique point A = Ao at which Fpp(q, Ao) matches the observable F(q). However, 
for cases (2.b) and (4. a) the situation is different but still there is a point A at 
which Fpxiq, Ao) can give a bitter result than that of perturbation theory. In all of 
these cases, we can estimate Ao from the properties of Fpx(q, A), or with the aid of 
higher order perturbative corrections, in a way that will improve on the perturbative 
predictions. On the other hand, in the remaining two cases (l.b) and (3. a) such an 
improvement does not occur. 

Having determined Ao as a function of q in a domain where perturbation theory is 
reasonably trustworthy, we use its average A e in FpT^q, A e ) to extrapolate to smaller 
values of q. This results in an estimate free of the unphysical divergences of the 
original perturbative approximation, and this is the central goal of our approach. 
This will be true even for the cases (l.b) and (3. a). In some cases, the effective 
average value A e is either equal to or not significantly different from Xo(q). This is 
always true if Ao(<?) is either a constant or a slowly varying function over a large 
domain of q. For instance, in our illustrative example A e assumes the value X ma x in 
case(2.b) and X min in case(4.a); in other words it coincides with A (g) for all q > A 2 . 
Also, it follows from the inequality: 

\r IR (q, A e )| < B¥£ e"^ 2 ) + M £±! for q ? A 2 , (2.39) 
q — A z tiA z A e q 

that the contribution of Tm(q, A e ) becomes negligible for q ^> A 2 . Hence, employing 
A e in the remote UV region does not spoil the standard perturbative results. 
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One final point to remark is that the method also applies in the reverse direction 
to the previous setting. For example, in the case when the infrared behaviour Fm(q) 
is well known then a similar representation to ( |2.1| ), but with q replaced by a variable 
with the dimensions of length squared r = 1/q, can be used in almost the same way 
as before except this time to deduce the corresponding analytic UV expression: 

F uv (q = l/r,\) = ^- [ ^1 Fi R (l/f)df. (2.40) 
2m J C IR r ~ r 

From this viewpoint, our method plays the role of a bridge between regions of small 
and large momenta, allowing us from a prior knowledge of either the UV or IR 
behaviour to extract information about the IR or UV properties respectively. 



3. Analysis of the one-loop QCD running coupling 

In this section, we demonstrate the implementation of our method in constructing a 
simple analytic model for describing the regular IR behaviour of the QCD effective 
coupling constant. For ease of calculation, we assume that the zeros of the exact 
function a(Q 2 ) of the running coupling are expected to occur only on the negative 
real axis Ke{Q 2 } < and/or infinity. This, together with the causality principle, 
allows us to express the reciprocal of ot(Q 2 ) in terms of the contour integral: 

= fc, j Cuv W) dk > (3 - 1} 

introduced in the previous section. When applying this representation to an n-loop 
perturbative expression aj^(Q 2 ), we take A = A e instead of oo. In this way, we 
obtain a combination of perturbative and non-perturbative contributions, which has 
the merit of: 

1. reinstating the causal analyticity structure in the right half plane Re{Q 2 }>0, 

2. preserving the standard UV behaviour and, 

3. improving the perturbative estimate in both the IR and low UV regions. 

In the following, we shall prove this claim within the one-loop approximation. For 
convenience, we express the reciprocal of aipp(Q 2 ) in terms of a dimensionless variable 
q = Q 2 /A 2 as: 

' £ Hi) , (3-2) 



aW(g) 47r 
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where a^(Q 2 /A 2 ) = ap^(Q 2 ) 2 . In the spirit of (|3.1|) , an analytically improved 
expression corresponding to the one-loop coupling constant can be defined as: 



a 



(1) (<?A) 



with 



X (1) 0?,A) 



1 

2tH 



„X(k-q) 



X=X e 



k-q aW(k) 



dk , 



(3.3) 



(3.4) 



Note that [a^ (k)]' 1 , being dependent on lng, has a cut along the negative real axis 
Re{k} < 0. We shall now proceed to evaluate the above integral for the case q ^ 0. 
Following the same argument of the preceding section from (|2.1|) to (|2.4j ), we can 
show that: 



X ^(q,\) = ^\n(q) + ^I Ck (q,\) 



ft 



47T 



where 



with 



W?,A) = J e (g,A) + ^" 



An 



Rw e -\(k+ q ) 



dk . 



J e (?,A) 



ee 



-\q 



2n 



e \ecos(e) Re 



k + q 

JXesin(e) \ n ^ e p *f 



e — qe 



-19 



dd. 



(3.5) 
(3.6) 
(3.7) 



A simple way to check that the contribution of J e (q, A) vanishes as e approaches zero 
is to consider the upper bound: 

p— X(q— e) 

|J«(ff,A)| < e — e(ln 2 (e) + 7r 2 ) 1/2 , 



q - e 



from which it follows that: 



lim J e (g,A) = 0. 



Consequently, in the limit as e — > we can write in a simple form: 



ic k (q,x) 



Rw e -X(k+q) 



(k + q) 

By inserting this into ( |3.5| ) and letting R vv 

A) 



An 



dk =E 1 (\q)-E 1 (\(q + R av )). 
oo, we obtain: 
\n(q) + E 1 (Xq) 



(3.8) 
(3.9) 

(3.10) 
(3.11) 



Having found the exact structure of x^(<?, A) for q ^ 0, let us now reconsider the 
contour integral ( |3.4|) once more to explore the continuity of x (<7> A) at the origin 
q = 0. Starting from the expression: 

„Xk 



X (1) (0,A) = ^-L 
Air 2m 



Cjj\ 



k 



ln(fc) dk . 



2 In this paper, we use the notation o^ n \q — Q 2 / A 2 ) = a^{Q 2 ) for the corresponding n-loop 
approximation. 



(3.12) 
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we arrive, after some work, at 



,">(0,AH§>im 



dk + 



ln(e) 



71 



3 Ae cos(8) 



cos(Ae sin(0)) dO 



1 

7T 



^Xe cos (9) 



sin(Aesin(0))0d0 



(3.13) 



Making use of 



and inserting the following series representations [|1(J 

e Aecos(e) cos ( Aesin (^)) cos (n0) for Ae < 1 



n=0 

oo 



7?! 



e A6co S (0) sin ( Aesin (#)) sin(n6) for Ae < 1 



n=l 



(3.14) 



into the two angular integrals in ( |3.13| ), we immediately obtain: 



* (1) (M) = ^Hm 



E 1 (\e) - E 1 (\R UV ) + \n(e) + J2 



4)n 



n=l 



n! n 



(Ae)" 



If we now make use of (|2.21|) and let R uv — > oo, we arrive at: 

X (1) (0,A) = ^ln(e-VA). 
This ensures the continuity condition at the origin: 

lim x (1 %A)=x {1) (0,A). 

Let us now find the range of A for which: 

a { >(q, A) 



X«(g,A) /3 Q [ ln(g) + Ei(Xq) 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



is continuous on the whole interval q > 0. This requires an investigation of the zeros 
of x (<7jA) in D^, where = {(g, A) : q > 0, A > 0}. From expression ( |3.11| ), 
we know that since Ei(Xq) > in D% and lng is negative only when q G (0, 1), the 
zeros of the function \ {Qj A) are expected to occur somewhere in the subinterval 
{(q, A) : < 5 < 1, A>0}c -D^. This domain can be reduced further by first 
considering the complementary exponential integral: 



Ein(x) 



Ml-e- 



-xt] 



dt = j2 



_^n-l x n 



—r , X > 



71=1 



n n! 



(3.19) 
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Figure 6: The zeros of A): the curve shown in the strip {(q, A) : < q < 1, A > e -7 } 

represents the roots of X (l-i^) in Dy.. This curve never touches the vertical q = 1 as 
X {li A) does not have a zero on this line. 



to rewrite ( |3.11|) as: 



X (1 \<1, A) = § [ ln(e"VA) + Ein{\ q) ] . (3.20) 

Then by observing the fact that Ein(Xq) > for all Xq > and that ln(e _7 /A) < 
only if A > e~ 7 , we deduce that the zeros of x (?> A) can exist only in the following 
interval {(q, A) : < q < 1, A > e -7 }. Hence, a^\q, A) is continuous for all q > if 
and only if A G (0, e -7 ). Another way to obtain this result is by considering the fact 
that a^(0, A) > a^(q, A) for all values of A that keep a^(q, A) finite and positive 
throughout the range q > 0. Then, from the positivity of: 

« (1) (0,A) = ^— i— -, (3.21) 
A) ln(e 7 /A) 

it follows that the allowed values of A are those confined to the interval: 

D eff = {A : < A < e" 7 w 0.561} . (3.22) 
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Using Maple VI for solving x (<Z> = numerically, we illustrate in Fig. (j^) the 
zeros of x^HQi i- e the singularities of a^(q, A), by a curve confined to the strip 
{(q, A) : < q < 1, A > e~ 7 }, showing that -D e // is the only region of A for which 
a^>(q, A) is continuous for all q > 0. 

4. The Estimation of A e 

In this section, we try to find an effective value A e G -D e // such that a^^g, A e ) 
provides a better estimate than that of the standard perturbative results, especially 
in the vicinity of the Landau-pole and the IR region. 

From the constraint ( |3.22|) , i.e. A e G D e ff, it follows that: 

< a (1) (g, A e ) < M(q) for all q > , (4.1) 

where the upper bound M(q) = a^^g, e -7 ). This implies that a^(q, A e ) < a^^g) = 
a^(q,\ 1) for all q > 1. However, since Ei(\ e q) in ( |3.18|) is negligibly small 
for q A e , we may consider a^(q, A e ) = a^(q) at sufficiently large values of q. 
For a given rif and g > 1, we find that the way in which a^i^q, A) changes with A, 
as shown in Fig.(|7|), is a typical example of case(l) illustrated in Fig.(f|). Hence, 
before setting up a criterion for determining A e , we need to find out whether the 
perturbative estimate a^^g), in the low UV region, overshoots or undershoots the 
true value a(g). This can be deduced by comparing a^(q) to the standard higher- 
loop corrections [17]]: 




Figure 7: The variation of a^(2,X) versus A, (a) for rif = 3 and (b) for < rif < 6. 
Here, a^(2,X) tends to a^(2) as A increases. In (b), aW(2, A) increases with rif for any 
fixed value of A. 
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Figure 8: The comparison between the 1-, 2- and 3-loop running coupling for rif = 3. 



and 



« (3) (<?) 



47T 



_/3j.lnL , A 2 ln 2 L-/?2 lnL + &#,-/3? 



where L = lng with g = Q 2 /A 2 , and: 



102 - 
2857 



38 



^ 71/, 

5033 325 , 
54 y 



(4.3) 

(4.4) 
(4.5) 



Fig.([8p illustrates the comparison in terms of the physical scale Q [GeV]. In this 
figure, we consider n/ = 3 and take A = 0.593 GeV, 0.495 GeV and 0.383 GeV for 
a PT(Q 2 )i a p%(Q 2 ) an d a i?%(Q 2 ) respectively. The A values are determined in the 



MS scheme from the experimentally measured r— lepton decay rate ]18[ given by: 



R T = 2.9087 (0.998 + 5^ n) ) = 3.492 , 



where 



5 {n) 
p 




4?(m 2 ; 



7T 



(4.6) 



(4.7) 



here n < 3 denotes the n th loop order approximation, ai = 1, = 5.20232, 03 = 
26.3659 and M T = 1.777 GeV. In Fig.(|8]), we used rif = 3 as the average number of 
active quarks, ignoring complications due to quark thresholds. This is reasonable in 
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n f 





1 


2 


3 


4 


5 


6 




2.204 


2.246 


2.293 


2.345 


2.402 


2.467 


2.539 



Table 1: Numerical values of qo{nf) 



the low energy interval Q < 3 GeV. For a more precise description of the evolution 
of a^l(Q 2 ) in a larger momentum interval, one should take into account the effects 
of quark thresholds. However, this does not change the overall picture illustrated 



in Fig.(|j). Having now observed that a^{Q 2 ) > or^(Q 2 ) > or^{Q 2 ), we would 
expect aS x >{q) to overshoot the true value of the coupling constant, at least for 
intermediate and sufficiently small values of q. This indicates that our problem is 
of the type classified earlier (l.a) and depicted in Fig.f^.l.a). Hence, we may 

ascertain that for every point q > 1 there exists a proper value \o(q) such that 
a^(q, X (q)) = «(<?), just as explained in case(l.a). One way to estimate Xo(q) is to 
first consider the curvature function of a^^g, A): 



K^(X,q) 



a 



(i) 

A A 



(?,A)| 



l + (a<£\q,\)y 



3/2 



(4.8) 



for a fixed q > 1, where 



a«(g, A) 
«aa(^ a ) 



^\q,X) = ^-[a^(q,X)} 



2 e 



-Xq 



d 2 



Arc 



X 



9A2 « (1) (9,A) = 4 1) (g,A)[^a( 1 )(g,A) ^ 



x~ q 



(4.9) 
(4.10) 



Then take Xo(q) to be the value of A at which K^(Xq, q) is a local maximum. For 
sufficiently small values of q > 1, K^'(X,q) has two local maxima, say at Ai and 
A 2 . In this case, we consider Ao to be the largest of these (say Ao = A2 > Ai) so 
that a^(q, A ) remains as close as possible to the perturbative results a^\q), which 
are still reasonable for such values of q. In Fig.(g) we plot the curvature function 
K^(X, q) versus A for some fixed values of q, illustrating in part (a) of the figure the 
difference between a^(q, A ) and a^(q) for q = 3. In part (b) of Fig.(||), we observe 
that the value of A decreases with increasing q. In general, if we plot A against 
q for any quark flavour nf we find that Xo(q) is a monotonically decreasing function 
as shown in Fig. (|10|) for nj = 3. Hence, we may express D e ff in terms of Xo(q) as: 



D eff = {A e : A (oo) < A e < A (g )} , 

7 . The values of go f° r different rif are all given in 



with q defined by X (qo) = e 
table(|l]). For sufficiently large values of q, our choice of Ao does not lead to a 
significant improvement on the perturbative estimate. In fact, since only the first 
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Figure 9: The variation of K ^(A; (?) versus A for nj = 3, (a) at q = 3 and (b) at q = 3, 4, 5. 
For a fixed q > 1, Ao(<?) denotes the value of A at which K^\X;q) is a local maximum. 
In (b), A (5) < A (4) < A (3). 
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Figure 10: The variation of Ao(g) versus g, for nj = 3. 



two decimal places of the coupling constant are of physical interest, we may take 
a^\q, Xo(q)) — a^(q) for all q > 120 and < nj < 6. On this basis, we shall focus 
our attention on the effective domain: 



Seff = { g ■ Qo<q< 120 } 



(4.11) 



over which the difference between a^'(q, Xo(q)) and a^\q) is appreciable. Thus, we 
expect to find A e in a narrower interval, D e ff = {A e : Aq(120) < A e < Ao(go)} C 



D e ff . Since, as illustrated in Fig. (|T0|), Xo(l) is a very slowly varying function over a 
large part of the interval S e ff, i.e. for q 6 (11, 120), we take the average of X (q) over 
S e ff to be the best estimated value for A e . Our numerical results for A e are listed in 
table(|2]) for the n/ values of practical interest. As seen from table(^), the values of 
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n f 





1 


2 


3 


4 


5 


6 




0.0379 


0.0390 


0.0402 


0.0416 


0.0429 


0.0444 


0.0461 



Table 2: Numerical values of A e (nj) 



A e do not change appreciably with n/. Hence, we shall fix A e to a central value of 
0.04 for all n f . 

Let us now denote a^(Q 2 /A 2 , A e ) by a (1 )(Q 2 ), where: 

aWfn 2 ^!^ (412I 

W J A, [ \n(Q 2 /A 2 ) + E 1 (X e Q 2 /A 2 ) ] ' 1 " ] 

In this expression, the exponential integral function E\(\ e Q 2 /A 2 ) plays an important 
role in removing the ghost pole at Q — A, preserving the correct UV behaviour (as 
it decays rapidly to zero for Q ^> A) and enforcing the running coupling to freeze 
in the low IR region. Therefore, it is essentially non-perturbative. At low energies, 
a^(Q 2 ) approximates well to the simple formula: 

5<1, « 2 '^^ ^ <Q^<^, (4,3) 

where n e = ln(e~ 7 /A e ) and to 2 = K e A 2 /A e . In particular, this yields a finite value at 
the origin, namely: 

Air 

a«(0) = -— , (4.14) 

P0«e 

which is totally independent of the characteristic mass scale A. Obviously, this is a 
key advantage of our approach as it agrees with the IR freezing phenomenological 
hypothesis jTT], [T2j . Now, let me remark from Ref . [12] that the freezing phenomenon 



is not incompatible with confinement as there is no evidence that confinement neces- 
sarily requires the coupling constant to become infinite in the IR region. For instance, 



Gribov's ideas |Cg, |13[ explicitly involve a freezing of the coupling constant at low 
momenta. This phenomenon is also present in perturbation theory but beyond the 
leading order and for certain numbers of quark flavours. 



5. Empirical Investigations of Our Model 

In this section, we compare our approach with other recognised methods including 
conventional and modified perturbation theory. Then, we test our model on phe- 
nomenologically estimated data for a fit-invariant characteristic integral depending 
solely on the IR behaviour of the coupling constant. 
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5.1 Threshold Matching 



In this subsection, we shall discuss the method we use for determining the values of 
A( n f> to be employed in our calculation. For sufficiently large values of Q 2 , our ex- 



pression for the coupling constant reduces to the familiar one- loop formula a^(Q 
This suggests that our must coincide with the usual QCD scale parameter A 
in momentum intervals corresponding to n/ > 5. Hence, we set our 5-flavour mass 



(n f ) 
PT 



parameter A^ 5 ^ equal to Ap^ at the energy scale of the Z-boson mass Mz = 91.19 
GeV. We extract the value A (5) = 0.135 GeV from a? T (M|) by means of the recent 
parametrisation of the hadronic decay width of the Z-boson : 



19.934 




(5.1) 



where n < 3 denotes the n th loop order under consideration, a\ = 1.045, ci2 = 0.94, 
a 3 = -15 and R z = 20.768 [271. To obtain the values of A^ and A {3 \ which may 



differ from Ap 4 ^ and Ap4 as they correspond to lower energy intervals, we use the 
familiar matching condition |p3fl : 

a {l \Q th - n f ) = a w (Q th ;n f - 1) , (5.2) 

at the energy thresholds Qth = 2mb and Qth = 2m c , where mi, = 4.21 GeV and 
m c = 1.35 GeV are the bottom- and charm-quark masses |ljj respectively. In this 
way, starting with A^ 5 ) = 0.135 GeV, we get A (4) = 0.188 GeV and A^ 3 ) = 0.229 
GeV. 

At light quark thresholds Qth < 2m s , where m s = 0.199 GeV is the strange- 
quark mass |TT 



.0) 



the applied matching condition ( |5.2j ) holds only approximately. 
This is also true for the analytic perturbation theory of Shirkov and Solovtsov ||. 
For simplicity, instead of introducing a nontrivial and more complicated matching 
procedure such as the one used in Ref . |pl| , we choose to take A^-* = A^ 3 -* for all 
n/ < 3. This choice is quite reasonable since it leads to: 



a (1) (Qt h ;n f ) 
aW(Q th ;n f -l) 

for all energy thresholds corresponding to «j < 3. 
ported by the good agreement with the value A^ ) ; 
from quenched lattice QCD in Ref. 



(5.3) 

Furthermore, it is also sup- 
= 0.238 (0.019) GeV obtained 
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5.2 Illustrative comparison 

We begin by considering a quick comparison between the predictions in our approach 
and those in perturbation theory at one-, two- and three-loop order. Fig.(|llJ) shows 
the behaviour of the 3-flavour running coupling a^(Q 2 ) together with ap^(Q 2 ), 
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a^Q 2 ) and af r {Q 2 ) in the low UV region as well as the Q-dependence of a^ l \Q 2 ) 
outside the perturbative domain. In this plot we use the same A values employed 
in Fig. (H) and take A^ 3 ) = 0.229 GeV for our formula a^(Q 2 ). As is seen from 
Fig.(|ll[b), our model improves significantly on the 1-loop perturbative predictions 
and agrees well with the higher- loop estimates throughout the range Q > 1.5 GeV. 
In Fig.([TT|.a), we observe that the evolution of a^^Q 2 ) slows down appreciably as 
Q enters the IR region, Q < A^ 3 \ and freezes rapidly to a finite value of 0.529. This 
provides some direct theoretical evidence for the freezing of the running coupling at 
low energies, an idea that has long been a popular and successful phenomenological 
hypothesis. Phenomenological studies, such as JTTJ and references therein, show 
that a running coupling that freezes at low energies can be useful in describing 
experimental data and some IR effects in QCD. In the low energy domain: 0.8 GeV < 
Q < 1.2 GeV , Badalian and Morgunov [p4]| extracted the value of the strong coupling 
constant a s (Q 2 ) from the fits to charmonium spectrum and fine structure splittings. 
They found that a s {Q 2 ) = 0.38 ± 0.03(exp.) ± 0.04(theory). This agrees reasonably 
well with our prediction, a^(Q 2 ) = 0.42 ± 0.03 in the same energy interval. 

Within the 3-loop approximation of the optimized perturbation theory, Mat- 
tingly and Stevenson [] 1| found that the 2-flavour QCD running coupling a s freezes 
below 0.3 GeV to a constant value of about 0.26 n. Although this value is quanti- 
tatively uncertain (while qualitatively unequivocal), it is not far off from our 1-loop 
prediction of 0.16 7r. 

A good theoretical approach which supports the IR freezing behaviour of our 
model follows from the work of Simonov and Badalian |3], |J. They studied the non- 
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n f 





1 


2 


3 


4 


5 


A^ /} GeV 
LO 


0.307 


0.283 


0.258 


0.230 


0.188 


0.135 


A^ } GeV 
NLO 


0.569 


0.553 


0.535 


0.509 


0.414 


0.274 



Table 3: The values of A„ / in LO and NLO 



perturbative contribution to the QCD running coupling on the basis of a new back- 
ground field formalism, exploiting non-perturbative background correlation functions 
as a dynamical input. In this framework, the one- and two-loop running coupling 
constants flU: 



a 



g>«« 



47T 



1 



A) ln[(<2 2 + M|)/A 2 



a 



Air 



(Q) = 7T 



Po ln[(Q 2 + M 2 )/A 



A ln[ln[(Q 2 + M 2 )/A 2 
(31 ln[(Q 2 + M 2 )/A^ 



(5.4) 



(5.5) 



were found to depend on a new parameter Mb (similar to A e in our model) known as 
the background mass. From a fit to the charmonium fine structure ||24|| , it was found 
that M B = 1.1 GeV. In table®, we list the values of for 

rtf < 5 in leading order 

(LO) and next-to-leading order (NLO). They are obtained as described in Ref.[B], 



using A^J, = 0.135 GeV in LO and A£f = 0.274 GeV in NLO, which is basically the 
same method that we used to obtain the values of our A*- n/ ). 

In Fig.(|l2|), we display the comparison between the 3-flavour coupling constant 
in our approach and those in the new background field formalism calculated to lead- 
ing and next-to-leading order. As shown in the figure, our low energy estimates lie 
approximately in the middle between the predictions of a^iQ 2 ) and a^iQ 2 ), in- 
dicating the self-consistency of our results. At the origin Q = 0, our results differ 
from that obtained by Shirkov and Solovtsov || by a small multiplicative factor 
k~ 1 = 0.38, i.e. aW(0) = 0.38 ain(0), on the other hand we find that the estimates 



.(5) 



of the background coupling constants (0) and a^ ; (0) are much closer to our pre- 



(2), 



dictions than to ain (0). Table(||) lists our results for a^(0) together with a4n (0), 



a^(0) and a^ ; (0) for < nf < 5. The numerical value of our coupling constant 
a^'(Q 2 ), as shown in table(^), is reasonably small which supports the notion of an 
expansion in powers of at low momentum transfers. Phenomenological verifi- 
cation of this fact would be of large practical value. Gribov theory of confinement 



.,(2), 



13] demonstrates how colour confinement can be achieved in a field theory of light 



fermions interacting with comparatively small effective coupling, a fact of potentially 
great impact for enlarging the domain of applicability of perturbative ideology to the 
physics of hadrons and their interactions [|9f . 
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(b) 



Q [GeV] 



Figure 12: Illustrative comparison of our formula a^(Q) with the LO and NLO expres- 
sions for the background coupling constants a^(Q) and a^(Q) for nf = 3 . 








1 


2 


3 


4 


5 


cv«(0) 


0.432 


0.460 


0.492 


0.529 


0.571 


0.620 


a^(0) 


1.142 


1.216 


1.300 


1.396 


1.508 


1.639 


4^(0) 


0.448 


0.448 


0.448 


0.446 


0.427 


0.391 


4 2) (o) 


0.713 


0.713 


0.714 


0.705 


0.576 


0.447 



Table 4: The values of a^(0) compared to those of ain (0), ct^(O) and o^(0) 



„( 2 ) ( 



In a number of cases of the QCD calculations it is necessary to estimate integrals 
of the form [ 
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a s (k 2 )i{k) dk. 



(5.6) 



where f(k) is a smooth function behaving like k p at k <C Q. In this integral, the 
interval of integration includes the IR region where the perturbative expression for 
the running coupling a s (k 2 ) is inapplicable. Hence, by introducing an IR matching 
scale Hi such that A <C hi ^ Q the contribution to integral ( |5.6| ) from the region 
k > fii can be calculated perturbatively. On the other hand, the portion of the 
integral below hi is expressed in Ref.p5| in terms of a non-perturbative parameter 
a p (Hi) as: 

a s (k 2 )k p dk 



(5.7) 



For hi — 2 GeV and p = 0, an excellent fit to experimental data yields: A(2 GeV) = 
a (2GeV)/7r = 0.18 ±0.03 § and A(2GeV) = 0.17 ±0.01 |§. These results agree 
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reasonably well with our estimate for A(2GeV), which is obtained from 
direct substitution of our 3-flavour formula ( f4.12|) for a s (k 2 ), giving: 



by 



Here \ P = 0.04 and A 



ln(fc 2 /A 2 ) +E 1 (X e P/A 2 
0.229 GeV. 



dk = 0.14. 



(5.f 



5.3 Gluon condensate 

In QCD instantons are the best studied non-perturbative effects, leading to the 



formation of the gluon condensate, an important physical quantity defined as [p8f| : 



K = (0\^ G%G« u \0) , 



71 



(5.9) 



where G a is the gluon field strength tensor and a s is the quark-gluon coupling 
constant. This can be related to the total vacuum energy density e t of QCD by 



means of the famous trace anomaly relation |27| : 







> 7T 



(5.10) 



where m q and denote the quark masses and spinor quark fields respectively. 
Sandwiching Q m between the QCD vacuum states in the chiral limit (i.e. m q = 0) 



and on account of the relation (Ol^lO) = Ae t [^8], P7| , one obtains: 
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K 



(5.11) 



In the dilute instanton-gas approximation with rif — 0, a direct relation between 
the gluon condensate K and the instanton density |28, |29| : 



p{r) 



2tt 



2tt 



cxp 



«(r- 2 ) 



a 



with r ^ 1/Q as the instanton scale size variable and b = 0.0015 



represented in the form [28 



K = 16 



5 dr ' 



(5.12) 



can be 



(5.13) 



where R c is a cut-off introduced by hand to avoid the uncontrollable IR divergences. 
Since our formula for the running coupling ( (4.12|) agrees reasonably well with the 
estimate of aprj, at short distances r < R c and is analytic for all r > 0, it is very 
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r [GeV -1 ] 




(b) 



r [GeV -1 ] 



Figure 13: The behaviour of the instanton density as expressed in terms of the running 
couplings (a) ctp^, and (b) 



tempting to use it in place of ctp^ to evaluate the integral in ( |5.13| ) in the limit 
R c — > oo. If we do this, we arrive at: 

f°° 1 1 



K = 0.25 x ll b 6 A / x b [\n(x- 2 )+E 1 (X e x 
'o 



exphyE^AeX- 2 )]^, (5.14) 



here A corresponds to rif = 0. A direct numerical integration of ( |5.14j ) yields: 

K = 4.8624 A 4 . (5.15) 

Using the value A = 0.229 GeV, which we computed in subsection (5.1) for small rif, 
we estimate the gluon condensate from ( |5.15|) as K = 0.013 GeV 4 . This is in good 
agreement with the value phenomenologically estimated from the QCD sum rules 
approach @, namely (0|^ G^GJJJO) = 0.012 GeV 4 . In addition, our estimate for 
the gluon condensate is also consistent with that calculated phenomenologically in 
Ref. 1301 and found to be K = 0.014 +°;°° 44 GeV 4 . Note that if we had considered the 
0.238 GeV, obtained from quenched lattice QCD in Ref. 



value A^ ) = 0.238 GeV, obtained from quenched lattice QCD in Ref. P%| , instead 
of our estimated value of A, we would have then found K to be 0.015 GeV 4 , which is 
still close enough to the two phenomenological estimates mentioned above. Hence, it 
may well follow that the instanton density p(r) is more likely to reveal more realistic 
and useful information when expressed in terms of our analytic running coupling 
In Fig.(^), we shed some light on the dependence of p(r) on the running couplings 



a 



(i) 



PT 



and which we shall denote by Ppt(^) and p{r) respectively. As seen in 
Fig.(|13|), at short distances, say r < 0.6 GeV" 1 , Ppt(^) and p(r) agree reasonably 
well whereas when we increase r to sufficiently larger values Ppt(^) rises very rapidly 
towards infinity while p(r) stabilizes to a fixed density. This reflects the practical 
usefulness of our approach to the strong coupling constant. 
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5.4 QCD (3— function and IR properties 

In quantum chromodynamics, the renormalisation group (3— function: 

P(a(Q 2 )) = Q 2 ^Jfp- ' (5-16) 



has the perturbative expansion [31 



* — \ / a \ 71+1 
P(a) = -aJ2 A. (^) , (5.17) 



n=0 



where /3 , (3\ and f3 2 are given by ( |1.1| ), (|4.4|) and ( |4.5|) respectively and fl32| : 

(3 3 = 29243.0 - 6946.30 n f + 405.089 n) + 1.49931 n) , (5.18) 

here the higher-order coefficients (3 n for n > 3 have not been calculated till now. 
Generally speaking, fl5.16j ) is the crucial equation that decides whether or not there 
is a so-called stable fixed point at the origin (i.e. whether f3(a) = for a = 0) of 
either the infrared or ultraviolet type, with implications for the absence or presence 
of asymptotic freedom in the gauge theory under consideration. The distinction 
between an infrared and an ultraviolet stable fixed point depends on whether the 
derivative (3' \a) > or f3\a) < at a = respectively. In both kinds of gauge 
theories (Abelian and non-Abelian) , there is a stable fixed point at the origin but it 
turns out that, while f3\a) > (stable IR fixed point) for an Abelian gauge theory, 
the reverse is true for a non-Abelian gauge theory , i.e. (3\a) < (stable UV fixed 
point). In general, it is extremely difficult to establish the existence of stable fixed 
points of a quantum field theory of (3{a) for q^O. In perturbative QCD, the one-, 
two-, three- and four-loop j3— function with quark flavours rif < 6 fail to exhibit any 
non-trivial zero that can be interpreted as a stable IR fixed point. Because of this, 
the perturbative QCD coupling constant «pt(Q 2 ) blows up as Q 2 — > 0. 

Our main task in this section is to demonstrate how the j3— function in our 
approach can provide useful information about the IR region that a truncated per- 
turbative series can not. As a cross check on our method, we shall calculate the 
stable IR fixed point for any number of quark flavour and show that it is consistent 
with our previous result (|4.14j ). Starting from ( [4.1 2|) with being renamed as a, 



our QCD (3— function in accordance with the definition (|5.16|) assumes the form: 



P(a) = -^a 2 [l-e- w ^], (5.19) 

where 

W{a) = e 0{Q) , (5.20) 
and the function 0(a) is defined by the transcendental relation: 

<P(a) + E 1 (e^ ) = p- + ln(A e ) . (5.21) 
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In the following domains of a: 



D l = {a: (b{a) >3} = {a:0<a< 2.021 flT 1 } 



(5.22) 



and 

D 2 = { a : <p(a) < -4 } = { a : 4.724 (3^ <a< 4.757 (3 Q l } 
a good approximate solution of Q5.21J ) can be found as: 

4:71 47T 

— + ln(Ae) - E 1 (X e exp (-5—)) if a E Di, 

p a Po a 



(5.23) 



a 



(5.24) 



if a E D 2 . 



We note that as a approaches zero, (f>(a) grows without limit, allowing (|5.19|) to 
reproduce the conventional one-loop QCD (3— function. It follows from ( |5.19| ) that 
for all rij satisfying Po( n f) > there exists a stable IR fixed point app identified 
by W(a-pp) = 0. Using the expansion of ( [2.21| ) in ( 5.21 ), we obtain, in the limit 
W(a) — > 0, the function: 

47T 

W{a) = - K e , (5.25) 



from which app is found to be: 



a FP 



foot 



47T 

P K e 



(5.26) 



This is completely consistent, as it should be, with our previous result (|4.14|) . 

For 6 < rif < 16, the higher-order terms of the QCD (3— function are known to 
permit the occurrence of IR fixed points. For example, for all values of rif such that 
(3o(rif) > and (3\{nf) < the two-loop (3— function f3p^(a) 3 possesses a non-trivial 
IR fixed point given by: 

(2 ) _ 4vr/5 



a 



FP 



01 



> 



(5.27) 



This may sound fine. However, the IR fixed points arising from the truncation of the 
perturbative series (|5.17| ) are likely to be spurious. For instance, at the candidate 



value for app, the first and second order terms in (3p^(a) are equal in magnitude 



x i.e. |/3o«fp /4vr| = app /(47r) 2 |), indicating the inefficiency of the two-loop 
approximation around the calculated IR fixed point ( |5.27|) . In fact, there is no way 
to prove that the IR fixed points extracted from perturbation theory alone are indeed 
the positive zeros of the true (3— function or anywhere near them. However, it is still 
worthwhile to test the perturbative expansion Q5.17 ) to 4-loop order against our 



3 The notations /3p^ (a) and a^p are used to denote the n-loop (3— function and its corresponding 
IR fixed point respectively. 



32 




-0.12 



(a) (b) 

Figure 14: The behaviour of the (5— function in our approach and in perturbation theory 
for rif = 6 (a) and rif = 9 (b). In (a), the curves of (3p^(a), (3p^(a) and f3p^(a) coincide 
and none of them achieves a perturbative IR fixed point. In (b), the locations of the fixed 
points are different but very close to each other. 



model for the (3— function. In Fig.(|n]), we demonstrate the full behaviour of our 
f3— function together with its perturbative counterpart for rif = 6 and 9. As seen in 
the figure, the (3— function in our approach (3(a) agrees very well with the correct 
perturbative predictions at and near the UV fixed point a = 0. As a leaves the 
region D, where D = {a : < a < 2.7/ /3 (nj) }, (3(a) starts to diverge significantly 
from ftp^a) and changes direction, turning back to zero in the limit a — > «fp as 
depicted in Fig.(|l4}a) for rif = 6. The occurrence of IR fixed points in /3p1p(a), 
(3prl(a) and (3p^(a) takes place only if rif £ [9, 16], [7, 16] and [8,16] respectively. In 
these rif— intervals, we find that the two-, three- and four-loop /3— functions behave 
qualitatively like /3(a). This is illustrated in Fig. (|T4"|.b) with (3prl(a) and (3p%(a) for 
rif = 9. It would have been interesting to see whether this qualitatively similar 
behaviour continues beyond the 4-loop order. Unfortunately, we know of no existing 
argument or calculation that definitively answers this question. In table([|), we list 
the values of the IR fixed points in our approach and those allowed in perturbation 
theory within the two-, three- and four- loop approximations for 6 < rif < 16. 

For rif G [8,16], the four-loop (3— function (3p^ has two positive roots. The 
smaller of these is an IR fixed point whereas the larger one is an UV fixed point. Thus, 
we shall denote the latter by au V FP- From table(f)]), we find that with increasing rif 

(2) (3) 

our IR fixed point «fp increases gradually whereas the perturbative roots a FP , a F p, 
a^p and a^ypp decrease to smaller values. However, for most values of rif we observe 
that the further a F p is from one root of some order, say Op P , the closer it is to another 
root of different order, say a^pp 1 ^. For example, for any rif e [8, 10] we find «fp to 
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n f 


7 


8 


9 


10 


11 


12 


13 


14 


15 


«FP 


0.751 


0.839 


0.951 


1.098 


1.297 


1.586 


2.039 


2.854 


4.757 


(2) 

a F p 


% 




5.236 


2.208 


1.234 


0.754 


0.468 


0.278 


0.143 


(3) 


2.457 


1.464 


1.028 


0.764 


0.579 


0.435 


0.317 


0.215 


0.123 


(4) 

a FP 


* 


1.550 


1.072 


0.815 


0.626 


0.470 


0.337 


0.224 


0.126 


(4) 
a UVFP 


* 


14.364 


12.090 


5.617 


3.294 


2.295 


1.781 


1.480 


1.286 



Table 5: The values of «fp compared to those of a FP , a F p and a F p. The * sign denotes 
the non-existence of a positive fixed point. 



be closer in magnitude to a F p and a FF than to a FP and a^y-pp- On the other hand, 
for any rif G [12, 14] our estimate for a FF is much closer to c^\pp than to any other 
value of the perturbative IR fixed points. Also, for rif = 11 we have a FF > ~ a FF . This 
analysis may not be conclusive however it provides us with some element of truth 
about our estimate for the possible magnitude of the true IR fixed point. 
It is noteworthy that (|5.26|) enables us to express A e in terms of a FF > as: 

A e = exp ( - 7 - ) . (5.28) 

\ PO "FP / 

This could have a direct beneficial effect on our approach if accurate phenomeno- 
logical estimates for the IR fixed points were available. Indeed, if this were the 
case Q5.28 ) would then allow us to determine the only free parameter A e in the the- 



ory straightforwardly. However, in the absence of this phenomenological data, the 
method we used in obtaining the value A e = 0.04 remains an adequate alternative 
as it results in a good approximation for the QCD f3— function with any number of 
quark flavours. 



6. Conclusions 

In this paper, we have developed a new approach for investigating physical observ- 
ables of the type F(Q 2 ) in regions that are inaccessible to perturbative methods of 
quantum field theory. In this approach, we have exploited the causality condition to 
reformulate F(Q 2 ) as a limit as A — > oo of a contour integral depending on A. In this 
way, we have shown that the perturbatively violated analyticity structure of physical 
observables can be reinstated in, at least, the right half of the complex Q 2 -plane by 
taking A = A e instead of oo, where A e > is the only free parameter of the the- 
ory. Explicit guidelines for finding an appropriate value for A e has been given with 
illustrative examples. We have also shown how our construction can play the role 
of a bridge between regions of small and large momenta, emphasising the fact that 
from a prior knowledge of either the UV or IR behaviour we can extract information 
about the IR or UV properties respectively. In fact, as our formalism incorporates 
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non-perturbative effects it extends the range of applicability of perturbation theory 
to cover the whole energy domain. 

We have demonstrated the implementation of our method in tackling the ghost- 
pole problem in QCD, giving a simple and new analytic one-loop expression for the 
strong coupling constant. As our approach incorporates an extra free parameter 
A e in addition to A, we have included a special method of determining A e from 
the UV behaviour of the A-dependent coupling constant a^(q,X). This method 
involves the calculation of the maximum curvature of a^(q, A) in the UV region 
q > 1 as explained and carried out numerically in section (4). In this framework, 
we have found the effective parameter A e to be 0.04 for 6i^(q,X e ). Although our 
prescription of obtaining A e via the maximum curvature procedure can be replaced 
by just tuning A e to allow the coupling fit with the experimental data available, 
we find our estimates for A e appropriate enough for our model, leading to a good 
agreement with the result a s (Q 2 ) = 0.38 ± 0.03(exp) ± 0.04(theory) extracted from 
the fits to charmonium spectrum and fine structure splittings |24]] at a low energy 
scale Q = 1.0 ± 0.2 GeV. A supplementary discussion of how to determine the QCD 
scale parameter A in our model has also been considered. 

A distinctive feature of our approach is that the running coupling freezes to a 
finite value at the origin Q 2 = 0, being consistent with a popular phenomeno logical 
hypothesis |TT], [12] namely the IR freezing phenomenon. This has also been supported 
by Gribov theory of confinement which demonstrates how colour confinement can be 
achieved in a field theory of light fermions interacting with comparatively small 
effective coupling, a fact of potentially great impact for enlarging the domain of 
validity of perturbative ideology to the physics of hadrons and their interactions. 

For illustration, we have carried out a comprehensive comparison between the 
predictions in our approach and those estimated by other theoretical methods. Be- 
sides, we have tested our model on a fit-invariant IR characteristic integral extracted 
from jet physics data ||, 25 1 over the energy interval Q < 2 GeV. This test shows 
a reasonable agreement with our prediction. For further applications, we have used 
our model for the running coupling to compute the gluon condensate and have shown 
that our result agrees very well with the value phenomenologically estimated from 
QCD sum rules P5| . Finally, we have calculated the (3— function corresponding to 
our QCD coupling constant and have shown that it behaves qualitatively like its per- 
turbative counterpart, when calculated beyond the leading order and with a number 
of quark flavours allowing for the occurrence of IR fixed points. 

In further studies, it would be interesting to apply our analytic approach to 
improve the standard expressions for the running masses in perturbative QCD. 
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